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Over 40 years after the calculation of the 1-loop QED corrections to the muon lifetime, new 
theoretical developments have made it possible to obtain an analytic expression for the complete 
2-loop QED contributions in the Fermi theory. The exact result for the effects of virtual and real 
photons, virtual electrons, muons and hadrons as well as e"'"e~ pair creation is 

Ar^^L = To {^y (ifl^ - ip^C(2) - ^C(3) + f C(4) + 53C(2) In 2 - (0.042 ± 0.002)) 

0^ . where To is the tree-level width. This eliminates the theoretical error in the extracted value of the 

CTN ' Fermi coupling constant, Gf, which was previously the source of the dominant uncertainty. The 

^•^ ' new value is 

J>^. Gf = (1.16637 ± 0.00001) x 10"^ GeV"^ 

<^ : 

^H The overall error has been roughly halved and is now entirely experimental. Several experiments are 

^H . planned for the next generation of muon lifetime measurements and these can proceed unhindered 

QQ ' by theoretical uncertainties. 

T-H ! 

CN ' I. INTRODUCTION 

>: 

l/^ . The three fundamental input parameters that enter into all calculations of electroweak physics are the electromag- 

^^ [ netic coupling constant, a, the Fermi coupling constant. Gp, and the mass of the Z^ boson, Mz- Their current best 
f^ ■ values, along with their absolute and relative errors are 

On" 

0> ; a == 1/(137.0359895 ±0.0000061) (0.045 ppm) 

^ ; Gf ^ (1.16639 ± 0.00002) x 10"^ GeV-^ (17 ppm) 

I ; Afz = 91.1867 ±0.0021 GeV (23 ppm) 

Qh. 
(D . 
r-| , In the mid-80's, just before the turn on of LEP, a GERN report concluded that the error on Mz would be ±50 MeV 

or 550 ppm and that "A factor of 2-3 improvement can be reached with a determined effort" . It was thus generally 

believed that the error on Mz represented the limiting factor in the precision with which theoretical predictions could 

be made. The situation has changed adiabatically, however, and the relative error on Mz now approaches that of Gp- 

C^ I Since LEP, as a machine, is not a radically new design, the lesson that we should take is that it is extremely difficult 

to predict the accuracy with which physical quantities will be measured, even in the relatively short term, and that 

one should constantly strive to reduce such errors to the minimum level consistent with the available technology. The 

possibility of precision physics at a muon collider serves to emphasize this point. 

With this in mind, and given the great cost and effort that was expended in reducing the error on Mz to its current 
value, it is reasonable to look again at Gp and see what is required to reduce its error to a level where it can never 
become an obstacle limiting the accuracy with which theoretical predictions can be made. 

Gf is extracted from the measured value of the muon lifetime, r^ = (2.19703 ± 0.00004) /is ffl and on the exper- 
imental side this is currently the source of the dominant error. New experiments are planned at the Brookhaven 
National Laboratory, the Paul Scherrer Institute and the Rutherford- Appleton Laboratory and it is likely that the 
uncertainty on Gp from this source will be reduced to somewhere in the range 0.5-1 ppm. 

Most of the work reported here appears in ref.s 



X 



II. THE FERMI COUPLING CONSTANT 

As given above the current relative error on the Fermi constant is dGp/Gp = 1-7 x 10^^. Of this 0.9 x lO^'^ is 
experimental and 1.5 x 10^^ is theoretical being an estimate of unknown 2-loop QED corrections. 
Gp is related to the measured niuon lifetime, r, by the formula 

- = r^-ro(i + Ag). (1) 

where 

Glml 

as calculated using the Fermi theory in which the weak interactions are described by a contact interaction. Aq 
encapsulates the higher order QED corrections and may written as a perturbation series in a^ = e^/(47r), the 
renormalized electromagnetic coupling constant. Thus 

oo 

A9-^A9« (3) 

i=0 

in which the index i gives the power of, a^ that appears in Ag^*^ . Note that Eq. (|l|) differs from the usual formula 0| 
in ways that begin to become important at the part-per-million level. It is known Mm that 



Ag(°) = -8a; - Ux^ In x + 8x^ - x'^ (4) 

A?^'^ = (v) (f - 3^(2)) + 0{arx\nx) (5) 

where x = ml/mi and C is the Riemann zeta function with C(2) = 7r^/6. That the Ag*^*^ remain finite in the limit 
rUe — > is a consequence of the Kinoshita-Lee-Nauenberg theorem [g[ whose discovery was largely prompted by this 
particular observation. 

Although the Fermi theory is not renormalizable, the Ag*^'^ can be shown [^ to be finite for all i. This remarkable 
feature follows from the fact that the V — A interaction is invariant under a Fierz rearrangement that interchanges 
the wavefunctions of the electron and the muon neutrino. Thus Fermi theory is equivalent to an effective theory in 
which the muon and electron occupy the same fermion current in the weak interaction lagrangian. After fermion 
mass renormalization is performed the divergences in the vector part of this current are independent of the fermion 
mass and hence cancel in exactly the way they would for QED. The lagrangian of Fermi theory is invariant under the 
transformations -ipe -^ 75'0e and rrig — > — Trig. The QED corrections to the axial vector of the part can thus be obtained 
from the those of the vector part by changing the sign of the electron mass and hence are finite as well. Moreover 
the two sets of corrections are equal in the limit rrie — > and, in that case, calculations need only be performed using 
the vector part of the Fermi interaction. This conclusion holds under any regularization prescription and avoids the 
complications associated with the use of 75 in dimensional regularization [ pO| . 

The foregoing discussion does not apply to the /3-decay of the neutron where the Fierz rearrangement generates 
scalar and pseudoscalar terms that bear no resemblance to QED and the radiative corrections are consequently not 
finite. 

III. THE 2-LOOP QED CORRECTIONS TO THE MUON LIFETIME 

The complete 2-loop QED corrections to the muon lifetime require the calculation of matrix element for the 
processes, fi^ —> i/f^e^De, ^J■^ —^ v^e~i>f^'^^ /i~ — > i'^e~Pe77 and /i^ -^ i'f^e~Vf^e^e~ with up to two virtual photons. 
All processes contain infrared (IR) divergences coming from either virtual photons, soft bremsstrahlung or both. 



The cancellation of IR divergences occurs between the various processes but this complication may be avoided by 
exploiting cutting relations and calculating the 2-loop corrections as imaginary parts of 4-loop diagrams, some of which 
are shown in Fig.s 2 and 3. In these Feynman diagrams thick lines represent a muon and the thin lines represent 
either the electron or the neutrinos all of which are taken to be massless. Since the external muon is on-shell any cut 
passing through a muon line will vanish and the only cuts contributing to the imaginary part are precisely the ones 
that generate the diagrams appearing in the calculation of muon decay. 

Recursion relations [Q obtained by integration-by-parts were first applied to reduce all dimensionally regularized 
integrals to a small set of relatively simple integrals. The well-behaved primitive integrals were then calculated by 
taking the external muon monrcntum, q, off mass shell to obtain expressions as power series in a; = —q'^/m'f^ and 
logarithms of x using well-established large mass expansion techniques Q . As the large mass expansion proceeds 
many terms, such as those that are topologically tadpoles, can be immediately discarded since they do not give rise 
to imaginary parts. Since the final result is required for x ^ 1 the complete series must be summed which can now 
be done in closed form in terms of polygamma functions and certain classes of multiple nested sums [|l3| . 

All diagrams were calculated in a general covariant gauge for the photon field and exact cancellation in the final 
result of the dependence on the gauge parameter was demonstrated. 

A. Hadronic Contributions 





(a) (b) (c) (d) 

FIG. 1. Hadronic contributions to muon decay after Fierz rearrangement of the contact interaction. 



Hadronic effects enter t^ at the 2-loop level through the diagrams shown in Fig.l. The shaded blob represents 
the hadronic vacuum polarization of the photon. The hadronic contribution can be calculated in the usual way 
using dispersion relations but, in contrast to other well-known situations for which such effects have been calculated 
|^,|5[, the momenta of the external fermions, to which the virtual photon is attached, is not fixed. Here the electron 
participates in the phase-space integration which complicates matters somewhat. 

Hadronic contributions are always afflicted to some degree by an uncertainty that arises from the experimental error 
on the measured cross-section (Thad = a{e'^e^ -^ hadrons). If this uncertainty turned out to be large there would be 
little point proceeding with the perturbative calculation of the 2-loop QED contributions to the muon lifetime. 



The shift induced in the inverse lifetime, F^, of the muon is given as a convolution integral 



AFhad = ^ / —R{mlz) AF(z) 

OTT J4p Z 



(6) 



over the hadronic spectrum, R{q^) = crhad/cpoint, and in which p = rn^/rn? = 1.61395... The convolution kernel, 
AF(z) is obtained exactly as an analytic function. When the integral is performed using actual hadronic data the 
result is 



AFhad = — ro 



(0.042 ±0.002) 



(7) 



which includes a rather conservative estimate of the hadronic uncertainty. Still the latter amounts to only 2 parts in 
10® and so is well under control. 

The integral (pi) can be used to obtain an expression for the contribution from diagrams where the hadronic vacuum 
polarization has been replaced by muon loop by setting 



2\ L 4 



i?Kz)= (^l + -j^l--. (8) 

which gives 

= -To ( — ] 0.0364333. (10) 

The result agrees with that subsequently obtained by perturbative methods. The effect of tau loops can be obtained 
in a similar way and, as expected on the basis of the decoupling theorem, is very small. 

B. Photonic Corrections 

Examples of photonic diagrams which when cut give rise to contributions to the muon lifetime at 0{a^) are shown 
in Fig.2. 





(a) (b) 





(c) (d) 

FIG. 2. Examples of diagrams whose cuts give contributions to /i~ ^ !/^e~i'e, ^~ ~+ i'fj,e~Uc.'y or /i~ -^ i^^e~!?e77. 

The result obtained for the complete set of photonic diagrams is 



(11) 



ar 



n 



2 



= ro — 3.55877 (12) 



where C(3) = 1.20206... and C(4) = vr^/gO. 

C. Electron-Loops and e^e Pair Creation 

Diagrams containing an electron loop whose cuts give contributions to muon decay are shown in Fig. 3. The result 
obtained for these diagrams is 





(b) 
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(c) (d) 

FIG. 3. Diagrams containing an electron loop whose cuts give contributions to muon decay, fi -^ v^e v^, /i 
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-r-'^' I W - S«^' - 5^'<=" 



= To ( — I 3.22034. 



(13) 
(14) 



The value given in Eq.([l4f) is consistent with a numerical study carried out by Luke et al. 16 1 in the context of 
semi-leptonic decays of heavy quarks. 

In order to obtain a UV finite answer the a diagrams in which the electron loop is replaced by the photon 2-point 
counterterm must be included and therefore a decision has to taken as to the renormalization scheme that is to be 
adopted. This will be discussed further in section IV. 



IV. THE RENORMALIZED ELECTROMAGNETIC COUPLING CONSTANT, an 

The use of dispersion relations to calculate the hadronic and muon loop contributions in the previous section 
naturally invokes a subtraction of the photon vacuum polarization at g^ = and is therefore equivalent to the on-shell 
renormalization scheme. In cases where there are two or more widely separated scales, such as rrie and ?7ip, use of the 
MS renormalization scheme is indicated since it automatically incorporates the large logarithms that arise into the 
value of the renormalized coupling constant, Ur, at tree level. 

It is therefore appropriate here to adopt the MS renormalization scheme. The hadronic contributions of section III. A 
that were obtained via dispersion relations must be corrected to convert them from the on-shell to MS renormalization 
scheme. As it turns out the contribution from muon loops is the same in both schemes when the 't Hooft mass is 
taken set to jjl — m^ as is appropriate here. It can be shown p^ that the MS renormalization scheme is implemented 
in a consistent manner by using the results of section III. A as they are given and setting 



ar 



Q.e{mi_, 






47r2 



2 

mf. 



2 " 



(15) 



where the logarithm of 0{a^) was first calculated by Jost and Luttinger |lg]. The substitution ( |15[ ) correctly resums 
logarithms of the form a" In"^^ (m^/rrig) for all n > and incorporates those of a'^ ln(m^/mg). Upon evaluation 



aeirrif,) = 1/135.90 = 0.0073582. 



V. CONCLUSIONS 

It has been over 40 years since the 1-loop QED corrections to the muon Ufetime were calculated. The 2-loop con- 
tributions have had to await the development of new theoretical techniques, as well substantial increases in computer 
speed and storage capacity, but are now available. 

The complete 2-loop QED contribution to the muon lifetime in the Fermi model may be encapsulated in the quantity 
Ag'^\ as defined in Eq.s(E]) and (0). After including the effects of virtual and real photons, virtual electrons, muons, 
taus and hadrons as well as e+e~ pair creation is found to be 

A,^^) ^ (^) ' (i||^ - ^C(2) - ^C(3) + f C(4) + 53C(2) In 2 - (0.042 ± 0.002)) (16) 

^^^"j (6.700 ±0.002) (17) 

The tiny effect of tan loops has been included in the numerical value given in eq. (|l 
This translates into a new value for the Fermi coupling constant of 

Gf = (1.16637 ± 0.00001) x 10"'^ GeV"^ (9 ppm) 

The error has been halved relative to its previous value and is now entirely experimental. 

New measurements of the muon lifetime are planned at the Brookhaven National Laboratory, the Paul Scherrer 
Institute and the Rutherford-Appleton Laboratory and it is therefore likely that the uncertainty on Gp from this 
source will be reduced to somewhere in the range 0.5-1 ppm. 

In that case the theoretical error should still be negligible but other issues such, as error on the muon mass, to^, 
and the upper limit on the muon neutrino mass, rrii, , need to be considered. 

Finally many of the results and techniques employed here can be readily taken over and applied to inclusive decays 
of the &-quark. 
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